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We report a series of electronic structure calculations for CrN using different exchange correlation 
potentials: PBE, LDA+[7, the Tran-Blaha modified Becke- Johnson, and hybrid functionals. In 
every case, our calculations show that the onset of magnetism in CrN should be accompanied 
by a gap opening. The experimentally found antiferromagnetic order always leads to an insulating 
behavior. Our results give further evidence that the Tran-Blaha functional is very useful for treating 
the electronic structure of correlated semiconductors allowing a parameter free description of the 
system. Hybrid functionals are also well capable of describing the electronic structure of CrN. 
The analysis of the system is complemented with our calculations of the thermopower that are in 
agreement with the experimental data. 
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I. INTRODUCTION 

Transition-metal nitrides (TMN) have drawn con- 
siderable attention due to their interesting physi- 
cal properties: mechanical strength, ultrahardness, 
corrosion resistance, high melting points, good elec- 
trical and thermal conductivities, etcX In addition, 
some TMN exhibit other properties such as super- 
conductivity (arising in VN (Ref. 0) and NbN 
(Ref. 3,4)), a transition from a superconductor to 
a Cooper pair insulator in TiN (Ref. 5), or a mag- 
netostructural phase transition (MPT) as a function 
of temperature in CrN (Ref. la). 

This MPT in CrN has attracted considerable at- 
tention and controversial results have been reported. 
Some studies have shown that this transition is 
linked to a discontinuity in the resistivity versus 
temperature curved— It has been characterized as 
a first-order MPT from a high-temperature para- 
magnetic (PM), cubic NaCl-type phase to a low- 
temperature antiferromagnetic (AFM) orthorhom- 
bic Prima one at a Neel temperature TV ~ 273- 
286 K£ The low-T AFM configuration is quite 
unique consisting of double ferromagnetic (FM) lay- 
ers stacked antiferromagnetically along the [110] di- 
rection (see AFM2 in the right panel of Fig. [T]). 
From neutron powder diffraction experiments, a 
magnetic moment of 2.36 /is per Cr atom was ob- 
tained by Corliss et al.^ while in Ref. a value 
of 3.17 hb was measured. This peculiar magnetic 
structure has been shown to be linked to the struc- 
tural distortion across TV i 12 ' 13 However, other works 
show the suppression of the phase transition due 



to epitaxial constraints when grown in films with 
the system remaining in the cubic phase at all 
temperatures J^r— 

The main focus of controversy has been the un- 
derstanding of CrN electrical behavior. On the one 
hand, the values of the resistivity at room temper- 
ature vary over several orders of magnitude (from 
mficm to kf2cm) in different works i£ ~ 10 i 14 ' 16 ~ ^£ On 
the other, the nature of the transition in terms of the 
electrical conduction behavior has had a wide range 
of interpretations in the literature: 

Metal to metal. In the work of Tsuchiya et al.^ 
polycrystalline CrNi +x films with 0< x< 0.2 have 
been grown by reactive sputtering. Nearly stoichio- 
metric films show an AFM first order phase transi- 
tion at TV ~ 260 K. Resistivity measurements for 
these samples indicated metallic behavior across TV 
with a steep decrease of the resistivity below the 
transition. The same behavior was found by Browne 
et al£^ in CrN powder. In the study of Inumaru et 
al.r^- CrN films were epitaxially grown on a-A^O^ 
(0001) and MgO (001) substrates by pulsed laser de- 
position. The films grown on C1-AI2O3 (0001) with 
its (111) planes parallel to the substrate show no 
structural transition whereas CrN(OOl) films grown 
on MgO (001) show a transition to an AFM state 
at about 260 K. Both films exhibit metallic In- 
dependence of the resistivity, but a drop in its values 
below the Neel temperature is observed for the sam- 
ples showing a MPT. 

Metal to semiconductor. Constantin et alJ grew 
CrN(OOl) thin films on MgO (001) using molecular 
beam epitaxy. The films show a MPT at 285 K, 



with activated semiconducting behavior above TV 
and a band gap of 0.07 eV. Below the transition tem- 
perature, the films are metallic. The polycrystalline 
samples of CrN prepared by Bhobe et al*£ show a 
MPT at 286 K. Above T N , CrN exhibits activated 
semiconducting behavior (with a gap of 0.07 eV) be- 
coming a disordered metal below TV- Evidences of 
strong electron-electron correlations were found, de- 
termining an on-site Coulomb energy (U) of 4.5 eV 
by resonant photoemission spectroscopy (PES). In 
Ref. [TH both single and polycrystalline CrN lay- 
ers have been grown by reactive sputtering depo- 
sition on both MgO and quartz substrates. Poly- 
crystalline CrN did show a disorder induced metal- 
insulator transition to an orthorhombic phase across 
T N ~ 280 K. 

Semiconductor to semiconductor. In the above 
mentioned work of Tsuchiya et ai.j2£ the non- 
stoichiometric films did not show a MPT. The resis- 
tivity decreases monotonically with T suggesting a 
semiconducting (hopping-like) behavior of the sam- 
ples. CrN powder samples synthesized by Herle et 
al*£ show a transition from a PM to an AFM phase 
at 280 K. Activated semiconducting behavior with a 
band gap of 0.09 eV was obtained from resistivity 
vs. T curves. Quintela et al£ studied polycrys- 
talline samples of Cri^V^N. Stoichiometric CrN 
shows a MPT at TV ~ 286 K and semiconducting 
behavior in the whole T-range. Above Tn, the resis- 
tivity vs. T curve shows activated semiconducting 
behavior with a band gap of 150 meV. Below the 
transition, the changes in the chemical bond linked 
to the AFM order, give as a result a non-activated 
semiconducting behavior. The thermoelectric power 
of stoichiometric CrN shows a linear T-dependence 
as expected for metals or highly degenerate semi- 
conductors. This supports the idea of an electronic 
state neither thermally activated nor fully itinerant 
in the AFM phase. Further evidence of this idea 
is given by the drastic reduction of the bulk modu- 
lus due to the bond softening and charge derealiza- 
tion into the shorter and less strong Cr-Cr bonds in 
the AFM phased By ultrahigh vacuum magnetron 
sputtering, Gall et a/j 14 i 16 ' 18 i 19 ' 24 grew single crys- 
talline CrN (001) films on a MgO(OOl) substrate. 
The absence of any discontinuity in the resistivity vs. 
T curves suggests that epitaxial constraints suppress 
the transition to an AFM phase with the film being 
cubic PM at all T. The resistivity decreases with in- 
creasing temperature indicating semiconducting be- 
havior over the whole T-range. From optical ab- 
sorption measurements, an optical band gap of 0.7 
eV was obtained ! 14 i 17 An indirect band gap of about 
0.19 eV was also estimated 24 - from optical measure- 
ments suggesting CrN is a Mott-Hubbard type insu- 
lator. In addition, the samples show semiconducting 



behavior in a variable range hopping regime below 
120 Ki&±2. 

As can be seen, there is a big discussion and 
disagreement between different experimental works. 
There are also some theoretical studies on CrN. 
Calculations within the local density approximation 
(LDA) have been performed on CrN by Filippetti 
et a/i 12 ' 13 For the distorted othorhombic structure, 
they found as ground state the experimentally ob- 
served AFM2 configuration. For the cubic struc- 
ture, the ground state is the so-called AFM1, sim- 
ilar to the AFM2, but with the spins antialigned 
every layer instead of every two layers (left panel 
of Fig. []}. This can be understood in terms of 
the magnetic stress: the AFM2 phase is under more 
stress with respect to the AFM1. This is due to the 
spin ordering asymmetry of the AFM2 phase on the 
(001) plane where each Cr has two spin-antiparallel 
and two spin-parallel nearest neighbors. Considering 
that bonds between the antifcrromagnetically cou- 
pled neighbors are under tensile stress, while bonds 
between ferromagnetically coupled atoms are under 
compressive stress, in the AFM2 magnetic config- 
uration, stress can be relieved by an orthorhombic 
distortion. In the AFM1 structure, it cannot be re- 
lieved because all the four-nearest neighbors are FM 
coupled. Hence, in order to relieve stress and re- 
duce the total energy, the system undergoes a shear 
distortion from cubic to orthorhombic. LDA calcula- 
tions show a metallic behavior for both AFM phases 
(weak metal since the density of states (DOS) at the 
Fermi level is small) . From studies of TM oxides and 
other correlated compounds it is known that LDA 
underestimates the Coulomb effects on the narrow 
d bands, predicting an unrealistic metallic state in 
systems that actually have a band gap. Hence, fur- 
ther electronic structure calculations in CrN within 
the LBA+U— method were performed by Herwad- 
kar et al££- If the on-site Coulomb repulsion (U) is 
considered, a band gap already opens for a value of 
U = 3 eV in the magnetic ground state, suggesting 
that CrN may be close to a charge-transfer insulator. 

To shed light on the existing controversies about 
CrN, we present the results of renewed ab initio cal- 
culations using different computational schemes to 
treat the moderately correlated d electrons of CrN 
and describe properly its electronic structure and 
transport properties. 



II. COMPUTATIONAL DETAILS 

Our electronic structure calculations were per- 
formed with the WIEN2k code ) 27 i 28 based on den- 
sity functional theory (DFT)22r— utilizing the aug- 
mented plane wave plus local orbitals method 
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FIG. 1: (Color online) Three dimensional plot showing 
the difference between the spin-up and spin-down (dif- 
ferent color) electron densities in the AFM1 (left panel) 
and AFM2 (right panel) phases. The typical d 3 (t2 9 ) elec- 
tronic configuration of each Cr 3+ cation and no magnetic 
polarization in neighboring ligands can be observed. 



(APW+lo)^ For the calculations of the transport 
properties we used the BoltzTraP code^ a semiclas- 
sical approach based on the Boltzmann transport 
theory. 

Most of the DFT-based calculations in solids have 
been done using two (semi)local approaches for the 
exchange correlation energy and potential: local 
density approximation (LDA)2ii or generalized gra- 
dient approximation (GGA). 34 Although they give 
useful total energies and structural parameters, both 
of them are known to underestimate the band gap 
of most semiconductors and insulators. 

There are several functionals that can provide 
more accurate values for band gaps. On the one 
hand, staying inside the true Kohn-Sham (KS) 
framework (multiplicative potential), the newly de- 
veloped semilocal potentials (such as the one de- 
veloped by Tran and Blaha based on a modifica- 
tion of the Becke- Johnson 3 - 5 - potential) can give good 
results! 36 ' 37 But on the other hand, nonmultiplica- 
tive potentials can also be used. These methods 
belong to the so-called generalized KS framework 3 - 
and most of them mix DFT and Hartree-Fock (HF) 
theories. The most used ones are the LDA+L7— ap- 
proach and hybrid methods. 39 

The LDA+17 scheme improves over GGA or LDA 
in the study of systems containing correlated elec- 
trons by introducing the on-site Coulomb repul- 
sion U applied to localized electrons (e.g 3d or 4/). 
We have performed calculations on CrN within the 
LDA+f/ method (using the fully localized version 
for the double-counting correction^ -) , taking U in a 
reasonable range for this type of 3c? electron systems 
(from 2 to 5 eV) comparable to the values obtained 
by PES experiments. The exchange parameter J 
was set to zero. 

In hybrid methods a certain amount of semilocal 
exchange is replaced by HF exchange with the cor- 
relation remaining purely semilocal. They usually 



lead to improved band structures in good agreement 
with experiments as shown for several semiconduc- 
tors and insulators or even strongly correlated elec- 
tron systems4ir— We have applied one of the most 
common hybrid functionals (PBEO)^^ to the study 
of the electronic structure of CrN. However, tradi- 
tional hybrid functionals are computationally very 
demanding due to the long-range nature of the HF 
exchange in solids. A way to reduce this problem 
of convergence is to consider only the short-range 
part of the HF exchange as proposed by Bylan- 
der and Kleinman4£ This scheme has been used by 
Heyd, Scuseria, and Ernzerho f 47 ' 48 who proposed a 
screened hybrid functional based on PBEO. We have 
also performed calculations using a screened hybrid 
functional as implemented in the WIEN2k code with 
the HF exchange screened by means of the Yukawa 
potential (so called, YS-PBE0)42 

The modified Becke- Johnson exchange potential 
(a local approximation to an atomic exact-exchange 
potential and a screening term) + LDA-correlation 
(from hereon TB-mBJLDA) allows the calculation 
of band gaps with an accuracy similar to the much 
more expensive GW or hybrid methods. 36 ' 37 ' 50 We 
have also studied the electronic structure of CrN 
by using the TB-mBJLDA potential which does not 
contain any system-dependent parameter. 

The calculations were done with R,MTK max = 7 
(the product of the smallest of the atomic sphere 
radii Rmt and the plane wave cutoff parameter 
K max ), which determines the size of the basis set. 
The chosen Rmt were 2.06 an 1.82 for Cr and N, re- 
spectively. The PBE, LDA+U, and TB-mBJLDA 
calculations were done with a 7 x 13 x 9 k-mesh 
for the integrations over the Brillouin zone. Since 
calculations with hybrid functionals are much more 
costly, a less dense k-mesh of 4 x 8 x 6 was used for 
the YS-PBEO and PBEO functionals. Convergence 
for the total energies presented is achieved with that 
k-mesh. The transport properties are very sensitive 
to the Brillouin zone sampling, therefore an even 
denser grid (19 x 37 x 27) was used to obtain con- 
vergence. 

As mentioned in Sec. U the transition at Tjy in 
CrN from the PM to the AFM phase is accompa- 
nied by a distortion of the rocksalt structure. The 
shear distortion consists of a contraction of the an- 
gle from the cubic (90°) to 88.23° (see Ref. OJ for a 
more detailed description). Calculations were done 
both, with and without distortion at the experimen- 
tal lattice constants (Ref. [a). In order to accomo- 
date the AFM2 phase, a unit cell comprising four 
formula units has been used. The calculations for 
the other magnetic phases were also performed with 
this unit cell. 
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III. RESULTS 

Since the compound is a moderately correlated 
TMN, a crude estimate of the possible electronic 
configuration of the material can be obtained from 
an ionic model. Taking the usual valency for N, 
the average valence for the octahedrally coordinated 
Cr is +3. In a high spin state (S = 3/2, mag- 
netic moment of 3 /its), the threefold-degenerate 
majority-spin channel t 2g is completely filled while 

the minority t% g remains empty with the e g states 
for both spin channels also completely empty above 
the Fermi level. This electronic configuration can 
be observed in Figfl] Depending on the balance be- 
tween crystal-field splitting (Acf) and Hund's rule 
coupling strenght {Jh), the first unoccupied bands 
could be e^ g (J H > A CF ) or t\ g (J H < A Ci r). Signif- 
icant contribution of the N p bands below the Fermi 
level is expected due to hybridization of these states 
with the Cr d bands. 

In order to analyze in more detail the electronic 
structure of the compound, we have performed cal- 
culations with various methods: PBE (which is of 
the GGA form), LDA+J7 with several values of 
U (2-5 eV), TB-mBJLDA, and hybrid functionals 
(screened and unscreened). Table U shows the rela- 
tive energies, band gaps, and magnetic moment on 
the Cr cation obtained for the AFM1, AFM2, FM 
and non-magnetic (NM) configurations. 

Looking at Table [H it can be seen that the or- 
thorhombic AFM2 phase is favored overall, inde- 
pendent on the functional. For the cubic phase, 
this magnetic arrangement is also the more favored 
for the LDA+U and hybrid functionals. The other 
phases are clearly less stable. For instance, with 
LDA+C AFM1 is disfavored by 35-41 meV/formula 
unit, FM by 146-209 meV/formula unit, and NM by 
1074-2456 meV/ formula unit. Hybrid functionals 
retain the relative energies between the AFM2 and 
AFM1 phases and make the FM phase even more 
unstable with respect to the ground state AFM2. It 
has to be remarked that for PBE within the cubic 
phase (as it occured with LDA in the calculations 
carried out by Filippetti et ah 12 ^ 13 ), the AFM1 phase 
is more favored than the AFM2 although by just 4 
meV/formula unit. An important point to note is 
that if the AFM1 is kept, no orthorhombic distortion 
would be predicted, independently of the functional 
used (i.e., the cubic AFM1 phase is always more sta- 
ble than the orthorhombic AFM1 one). This is a 
consequence of the magnetostriction/stress relief is- 
sues described above and already found by Filipetti 
et al. using LDAH22. 

Concerning the band gap values, both the NM and 
FM configurations remain metallic for all the func- 



tionals except PBEO with a = 0.25, which gives a 
band gap around 0.45 eV for FM. In LDA+U calcu- 
lations with U = 3 eV (both for cubic -undistorted- 
and orthorhombic -distorted- structures), a gap of 
0.3 eV already opens for the AFM2 configuration 
while the AFM1 remains metallic for this particu- 
lar value of U (a gap of about 0.2 eV opens for the 
AFM1 ordering at U = 4 eV). If the value of U is 
increased up to 5 eV, the value of the band gap is 0.9 
eV for the AFM2 configuration and about 0.5 eV for 
AFM1. Hence, from LDA+U calculations it can be 
concluded that the opening of a band gap is easier 
for the AFM2 (ground state) phase. The same con- 
clusions can be drawn from the TB-mBJLDA cal- 
culations: both AFM configurations are insulating 
with values of the band gap similar to those obtained 
within LDA+U for U = 4 eV (0.79 and 0.22 eV for 
the AFM2 and AFM1 distorted structure, respec- 
tively). For hybrid functionals, the values of the 
band gap are very dependent on the amount of HF 
exchange (a) and on whether it is screened or not. 
For a = 0.25, PBEO gives a band gap for the AFM2 
distorted phase of about 2.1 eV and of 1.5 eV for 
AFM1, whereas the screened YS-PBE0 gives much 
lower values of the band gap (1.45 eV for AFM2 and 
0.84 eV for AFM1). Furthermore, the reduction of 
a to 0.1, leads to a huge reduction of the gap val- 
ues, making the AFM1 phase metallic and reducing 
significantly the AFM2 band gap. 

Concerning the magnetic moments of Cr ions, the 
values for LDA+C/ are in between the reported ex- 
perimental values of 2.36 fisr^ and 3.17 [Ib^ and in- 
crease with U. They change from 2.51 to 2.69 fis for 
the AFM2 configuration and from 2.43 to 2.62 hb for 
the AFM1 when varying U from 2 to 5 eV. Slightly 
larger values are obtained for the FM phase (from 
2.61 to 2.76 hb), reaching values close to the atomic 
limit for a d 3 cation in the high-spin state. This in- 
crease is due to the reduction of the p-d hybridization 
caused by U. For TB-mBJLDA, the magnetic mo- 
ments are close to the values obtained from LDA+U 
with low values of U = 2-3 eV. Hybrid functionals 
(both screened and unscreened) give magnetic mo- 
ments which are slightly larger (smaller) than TB- 
mBJLDA for a = 0.25 {a = 0.1). 

Below, a more detailed discussion of the electronic 
structure is given for each functional separately. 



A. LDA+C7 electronic structure analysis 

The evolution of the band structure with the U 
value for both the cubic AFM1 and orthorhombic 
AFM2 phases can be seen in Fig. [3] We use an 
orthorhombic Brillouin zone in all cases and the fol- 
lowing k-path: M(ir/a, 0, 0)-r(0, 0, 0)-X(0, n/b, 
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TABLE I: Relative energy (in meV/formula unit), fundamental band gap (in eV), and Cr atomic moment (in hb) 
of CrN for the different methods used in this work. It was not possible to converge the calculations with the hybrid 
functionals for the metallic NM phase. 
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Q)-W(it/a, 3tt/2&, 0)-L(0, n/b, n/c)-T-Z(0, 0, tt/c). 
In Fig. [3J the density of states (DOS) for the or- 
thorhombic AFM2 phase (for [7 = 4 eV) is also 
shown. 

For both magnetic configurations, the rough ionic 
electronic structure described above for the majority 
spin state of Cr atoms can be easily distinguished: 
the t 2g manifold is filled (in the range between —2 
and eV below the Fermi energy) , while only a small 
amount of (bonding) e g levels (in the range between 
—6 and —3 eV and hybridizing with the N p or- 
bitals) are occupied. As can be seen in Fig. [3l the 
band gap character is clearly d-d (i| g -ej) although 
some rather small contribution of N p states is also 
observed. Hence, CrN can be better described as 
a Mott-Hubbard insulator than as a charge-transfer 
insulator. 51 

For the AFM2 configuration (upper panels of Fig. 



[2]), it can be seen how for the lowest U value, the 
band structure remains metallic with a gap opening 
already for U = 3 eV. The unoccupied bands move 
up gradually as U increases. This is particularly 
noticeable for the first unoccupied band at T for U = 
2 eV and also the one that forms the small electron 
pocket. Also, as U increases, the bands right below 
the Fermi level move down in energy. For larger U, 
the Cr d character near M is reduced since the Cr 
t 2g states arc shifted down in energy. The splitting 
between the first two bands right below the Fermi 
energy decreases at V with increasing U, increases at 
M and W causing in principle a displacement of the 
top of the valence band from W for lower U values 
to M for higher ones (for U = 4 eV, almost the 
same eigenvalue is obtained at M and W). Some 
features of the band structure are U-independcnt: 
the conduction band minimum is always at T and 
the existence of degenerate bands at L which is not 
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FIG. 2: (Color online) Band structures with band character plot (Cr d highlighted) of the orthorhombic AFM2 phase 
(upper panels) and cubic AFM1 (lower panels) of CrN obtained from the LDA+C/ method with (a) U = 2 eV, (b) 
U — 3 eV, (c) U = 4 eV, and (d) U = 5 eV. The colors indicate the irreducible representation of the eigenvalues. 



accidental but symmetry required. 

Another observation that can be made by look- 
ing at the changes in the DOS, is a decrease of the 
distance between the main peaks of the occupied Cr 
tig and N p bands. For U = 2 eV, the distance is 
approximately 4 eV, while for U = 5 eV it is reduced 
to about 2.5 eV. 

Previous band structure calculations^ reported 
that when the gap first opens at low U, the top of the 
valence band (VBM) is at W moving to M as U in- 
creases. Figure @] shows the LDA+C/ band structure 
along the path M (w/a, 0, 0)-W (w/a, 3ir/2b, 0)- fc.l 
(ir/a, 3tt/26, tt/c)-W (tt/o, tt/26, 0). Our calcula- 
tions show that when the gap opens (at U ~ 3 eV), 
the VBM is not at W but displaced to a non-high 
symmetry point in the direction from fc.l to W' , and 
then moving to M at U = 5 eV. 

For the cubic AFM1 phase, a small gap of 0.20 eV 
is opened for U = 4 eV becoming 0.56 eV for U = 5 
eV. We can see this evolution in the lower panels 
of Fig. [21 In addition, numerous degenerate bands 
can be seen at M and W (with these degeneracies 
clearly lifted in the orthorhombic AFM2 phase). 



B. TB-mBJLDA electronic structure analysis 

As mentioned above, the band gap values obtained 
for TB-mBJLDA are rather similar to those ob- 
tained from LDA+C/ with U = 4 eV (Bhobe et al^. 
determined an experimental value of U of about 4.5 
eV by PES). 

For the cubic AFM1 configuration, the band 
structure plot is shown in Fig. [5] (right panel). Be- 
low the Fermi level it is fairly different from the one 
obtained within LDA+C/ (with low U). In addi- 
tion, the first band above the Fermi level at T in 
LDA+C/ calculations lies lower in energy than with 
TB-mBJLDA, but it is not just a rigid band shift, 
some differences are noticeable particularly around 
the r point. Actually, the bottom of the conduc- 
tion band in TB-mBJLDA is no longer at T but dis- 
placed to Z. In addition, the band gap character 
changes from LDA+C/ results being in this case a 

tig-tig band g a P- 

For the orthorhombic AFM2 phase, the splitting 
at r between the first two bands above the Fermi 
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FIG. 3: (Color online) Density of states of CrN 
(orthorhombic AFM2 phase) calculated with different 
methods. The Fermi energy is set at E = eV. The 
positive and negative values on the y-axis are for spin- 
up and spin-down, respectively. 



level is reduced in TB-mBJLDA (see left panel of 
Fig. [5]). The band structure obtained below the 
Fermi energy is very similar to the one obtained from 
LDA+C7 with U = 2 eV. In addition, as for LDA+U 
calculations with the values of U 3 to 4 eV, the VBM 
is not at W but displaced to a non-high symmmetry 
point in the direction k.l — W'. In Fig. [3] we can 
compare the DOS calculated from LDA+t7 and TB- 
mBJLDA. We can see that the N p character in the 
Cr t\ g band (from —2.5 to eV) is smaller in TB- 
mBJLDA than in LDA+U (and the opposite is true 
for the Cr character in the N p bands below -3 eV) . 
This effect has to do with the different hybridizations 
obtained with the two schemes (enhanced with the 
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FIG. 4: (Color online) Band structures blown-up around 
Ef for the orthorhombic AFM2 phase within LDA+U 
[(a) U = 2 eV, (b) U = 3 eV, (c) U = 4 eV, and (d) 
U = 5 eV] in a special k-path showing the displacement 
of the VBM with the U value. The colors indicate the 
irreducible representation of the eigenvalues. 



LDA+U method). In addition, the separation of N p 
and Cr t\ g bands is larger in TB-mBJLDA, although 
the total valence band width is almost the same (~ 
7.5 eV) 

The optical transmission and reflection spectra of 
CrN are shown in Ref. [It]- A band-gap-like op- 
tical transmission feature is found at 0.7 eV sug- 
gesting that CrN is a semiconductor with an indi- 
rect gap of approximately 0.7 eV and a T- valley gap 
greater than 0.7 eV. The values of the band gap 
(indirect) obtained from LDA+C7 (U = 4 eV) and 
TB-mBJLDA calculations for the AFM2 phase are 
very close to this 0.7 eV value (0.68 eV and 0.79 eV, 
respectively). A gap at T of about 1 eV (> 0.7 eV) 
is also found in our calculations. 

Thus, TB-mBJLDA is able to reproduce some of 
the effects on the electronic structure caused by the 
on-site Coulomb repulsion between Cr d electrons 
without using a tunable parameter and being con- 
sistent with the results found for bulk CrN stoichio- 
metric samples. 



C. Hybrid functionals electronic structure 
analysis 

As can be seen from Table U the results with hy- 
brid functionals strongly depend on the amount of 
HF exchange a. If we take as a reference the val- 
ues obtained within TB-mBJLDA (similar to those 
from LDA+f/ with U = 4 eV), using a = 0.1 gives 
closer results than using a = 0.25. The band struc- 
ture for the cubic AFM1 and orthorhombic AFM2 
phase are shown in Fig. [6]for the screened YS-PBE0 
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FIG. 5: (Color online) Band structures with TB- 
mBJLDA for the orthorhombic AFM2 (left panel) and 
cubic AFM1 (right panel) phase. The colors indicate the 
irreducible representation of the eigenvalues. 

functional with a = 0.1. 

For the ground state orthorhombic AFM2 phase, 
there is a reduction of the band gap with respect 
to LDA+J7 or TB-mBJLDA (as a consequence of 
the reduction of a) to 0.2 eV. The band structures 
retain the same features that we discussed above for 
LDA+J7 or TB-mBJLDA. In particular, the VBM is 
also not at W, but displaced to a non-high symmetry 
point (see Fig. 2] for LDA+U). Above the Fermi 
level, the splittings between the bands show more 
similarities with those obtained from LDA+C/ with 
U = 4 eV. Figure El shows that the shape of the DOS 
is very similar to the TB-mBJLDA DOS, but the 
total valence band width is increased to about 8.2 
eV. The same thing can be observed for the AFM1 
phase: the band structure below the Fermi energy is 
more similar to the one obtained with TB-mBJLDA 
and above to that obtained with LDA+t/. 

Hence, hybrid functionals can also describe accu- 
rately the electronic structure of CrN, but by tuning 
the amount of HF exact exchange. Although results 
and tunability with hybrid functionals are compara- 
ble with LDA+C/, this last method allows an easier 
physical description of the tuning parameter U and 
is computationally less expensive. 



IV. DISCUSSION 

As can be seen from the values of the band gap in 
Table HI the tendency towards a band gap opening 
is strongest for the AFM2 ground state phase (for 
LDA+U with U as low as 3 eV, TB-mBJLDA, and 
hybrid functionals). Although more difficult, it is 
also possible to open a gap in the other AFM order- 
ing (AFM1). 

As mentioned in Sec. HI the intrinsic mechanism 



FIG. 6: (Color online) Band structures with YS-PBEO 
(a = 0.1) for the orthorhombic AFM2 (left panel, in- 
sulating) and cubic AFM1 phase (right panel, metallic). 
The colors indicate the irreducible representation of the 
eigenvalues. 



of electronic transport in CrN is a controversial is- 
sue. Most of the discrepancies in experimental works 
appear among those where the properties of thin 
films (where by epitaxial constraints even the phase 
transition can be suppressed) have been measured. 
Hence, we will focus in the couple of works where 
bulk polycrystalline samples have been studied to 
reconcile our results with previous data^^ 

In the work by Quintela et al£ the resistivity of 
Cri_ x V x N has been measured. Negative tempera- 
ture coefficient of resistivity (TCR = dp/dT), typ- 
ical for semiconducting behavior, was observed for 
samples with V doping x < 0.1. For x > 0.2 the 
system shows a positive TCR pointing to an itiner- 
ant electron behavior. Stoichiometric CrN exhibits 
TCR < (nonitinerant) over the whole T-range al- 
though only for T > Tjy the behavior can be fit- 
ted to that of an activated semiconductor with a 
value of the activation energy of 75 meV. In the low- 
temperature AFM orthorhombic phase, the mech- 
anism of electronic transport is not fully under- 
stood: neither a simple thermal activation nor a 
variable range hopping behavior describe the resis- 
tivity under Tjv- In a system with a charge gap, 
given the exponential dependence of p(T) with T, 
d(ln a)/d(ln T) should diverge as T — > 0. However, 
for CrN, it goes to zero as T — s> 0. Hence, despite 
the negative TCR, CrN cannot be classified as a typ- 
ical thermally activated semiconductor. The differ- 
ent semiconducting behavior in the low-T phase was 
ascribed to the changes in chemical bond linked to 
the AFM order. Deriving the same result for stoi- 
chiometric CrN, Bhobe et alM> concluded that the 
low temperature phase of CrN is one of the rare ex- 
amples of AFM-itinerant (metallic) systems in spite 
of negative TCR. Back to Ref. 0, for x = 0.1 in 



the Cri-^V^N series, there is a sudden change in 
the sign of TCR at the Neel temperature: TCR > 
(itinerant) for the high-T PM phase and TCR < 
(nonitinerant) for the low-T AFM phase. Hence, it 
seems that for T < Tjv, the opening of a charge gap 
at the magnetic ordering temperature increases the 
activation energy, recovering semiconducting like be- 
havior below T/v- 

The fact that the opening of a charge gap is intrin- 
sic to the ground-state AFM ordering is consistent 
with our calculations. We show that the onset of 
magnetism should be accompanied of a gap open- 
ing. The behavior of the transport properties in the 
PM phase is beyond the scope of this work, but our 
calculations can conclude that the AFM phase is al- 
ways insulating and, if it is not, then it has to do 
with issues caused by nanostructuring such as non- 
stoichiometries, strain, surface effects, etc. In a sys- 
tem situated close to a metal-insulator transition, 
small changes can induce a transition from localized 
to itinerant behavior, but this is not reachable in 
bulk-stoichiometric CrN in the magnetic phase, as 
our calculations show. 

Further evidence of the existence of an unconven- 
tional electronic state in the AFM phase of CrN, 
not thermally activated neither fully itinerant, was 
observed in Ref. comparing the behavior of the 
electrical resistivity with the thermoelectric power 
measurements. For stoichiometric CrN, the thermo- 
electric power shows linear temperature dependence 
(typical for metals or highly degenerate semiconduc- 
tors). In addition, no discontinuity is observed at 
x = 0.2, where TCR changed from negative to posi- 
tive at low temperatures. 

In order to check if the experimental behavior of 
the thermoelectric power can be reproduced from 
our calculations, we have obtained its dependence 
with temperature in the magnetic ground state 
(AFM2) for LDA+/7 (U = 4 eV) and TB-mBJLDA 
calculations. This has been carried out from our 
band structure calculations within a semiclassical 
approach based on the Boltzmann transport theory 
through the BoltzTraP code^ 

For the sake of comparison, we present the calcu- 
lations together with the experimental values of the 
Seebeck coefficient taken from Ref. 1 (see Fig. ED. 
The results for the AFM2 phase fit the experimen- 
tal values rather well, particularly in TB-mBJLDA, 
both in order of magnitude of the thermopower and 
also on the observed evolution with temperature 
(dS/dT < 0) (although the behavior found from the 
calculations is not as linear as obtained experimen- 
tally). 

Hence, our thermopower calculations confirm that 
stoichiometric CrN in the magnetic phase is semi- 
conducting, but close to a transition from localized 
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FIG. 7: (Color online) Experimental 9 and calculated 
temperature dependence of the thermopower. The re- 
sults for both TB-mBJLDA and LDA+U (U = 4 eV) 
are shown. 



to itinerant electronic behavior. 



V. SUMMARY AND CONCLUSIONS 

We have performed electronic structure calcu- 
lations for CrN using a diverse set of exchange- 
correlation potentials: PBE, LDA+C/, the semilocal 
functional developed by Tran and Blaha based on a 
modification of the Becke Johnson potential,— and 
hybrid functionals both in screened (YS-PBEO) and 
unscreened (PBEO) modes. 

In any case, our calculations show that the onset 
of magnetism in CrN should be accompanied of a gap 
opening. We can conclude that the AFM phase is 
always insulating in agreement with experiments in 
stoichiometric bulk samples. The different metallic 
behavior found in various thin films has to do with 
nanostructuring that could include vacancies, strain, 
surface effects, etc. In a system situated so close 
to the itinerant boundary, small changes can yield 
mctallicity and even the suppression of the magnetic 
transition (which we have seen from our calculations 
it is strongly dependent on the structural details). 
The behavior obtained for the thermoelectric power 
confirms that stoichiometric CrN in the magnetic 
phase is semiconducting, but close to the transition 
from localized to itinerant behavior. The solution 
we found as ground state can describe accurately 
the transport properties of the system. 

In addition, our results give further evidence that 
the TB-mBJLDA functional is very useful for treat- 
ing the electronic structure of correlated semicon- 
ductors giving almost the same band gaps as ob- 
tained from LDA+C7 with U = 4 eV (close to the 
one determined experimentally) for different mag- 
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netic configurations allowing a parameter free de- 
scription of the system. Hybrid functionals are also 
well capable of describing the electronic structure 
of CrN by tuning the amount of HF exact exchange. 
However, results and tunability are comparable with 
LDA+C/, that allows an easier physical description 
of the tuning parameter U. 
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